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Abstract 

The differential cross section for the inclusive production of isolated prompt photons 
has been measured as a function of the photon transverse energy Ej in pp collisions 
at v^s = 7 TeV using data recorded by the CMS detector at the LHC. The data sample 
corresponds to an integrated luminosity of 2.9 pb^^. Photons are required to have a 
pseudorapidity \rj^\ < 1.45 and Ej > 21 GeV, covering the kinematic region 0.006 < 
xj < 0.086. The measured cross section is found to be in agreement with next-to- 
leading-order perturbative QCD calculations. 

Submitted to Physical Review Letters 




*See Appendix A for the list of collaboration members 



1 



The measurement of isolated prompt photon production in proton-proton (pp) coUisions pro- 
vides a test of perturbative quantum chromod3mamics (pQCD) and the possibility to constrain 
the parton distribution functions (PDF) of the proton. Such a measurement complements deep- 
inelastic scattering, Drell-Yan pair production, and jet production measurements [1-3|. At the 
Large Hadron Collider (LHC) IH, a significant increase of centre-of-mass energy with respect 
to previous collider experiments |'5'-9| allows the exploration of new kinematic regions in the 
hard-interaction processes in hadron-hadron collisions |10|. Isolated prompt photon produc- 
tion also represents a background to searches for new phenomena involving photons in the 
final state. 

In high-energy collisions, single prompt photons are produced directly in qg Compton scat- 
tering and qq annihilation, and in the fragmentation of partons with large transverse momen- 
tum. Photons are also produced in the decay of hadrons, which can mimic prompt production. 
This background, mostly from the decays of energetic and rj mesons, can be reduced by 
imposing isolation criteria on the photon candidates. 

This Letter presents a measurement of the differential production cross section of isolated 
prompt photons as a function of the photon transverse energy E j in pp collisions at a/s = 
7 TeV. The analyzed data sample corresponds to 2.9 ± 0.3 pb^^ of integrated luminosity, as 
recorded by the CMS detector at the LHC [llj. Isolated prompt photons with a pseudorapidity 
< 1.45 and > 21 GeV are studied. Here, //^ = - ln[tan(0/2)] and £^ = sin(0), 
where E"^ is the photon energy and 6 is the polar angle of the photon momentum measured 
with respect to the counterclockwise beam direction. This measurement exploits the differ- 
ence between the electromagnetic shower profiles of prompt photons and of photon pairs from 
neutral-meson decays. 

Photons are detected in the lead tungstate (PbW04) crystal electromagnetic calorimeter (ECAL), 
covering \tj\ < 3.0, comprising barrel and end cap sections. The analysis presented in this Let- 
ter is restricted to the barrel section, which covers \ri\ < 1.479. Light produced in the crystals 
is read out by avalanche photodiodes (APD) in the ECAL barrel. The ECAL barrel granu- 
larity is Afj X A(p = 0.0174 x 0.0174, where (p is the azimuthal angle measured with respect 
to the beam direction. The ECAL has an ultimate energy resolution better than 0.5% for un- 
converted photons with transverse energies above 100 GeV ||T2|. In 2010 collision data, for 
Ej > 20 GeV, this resolution is already better than 1% in the barrel flS]. Surrounding the ECAL 
there is a brass and scintillator sampling hadron calorimeter (HCAL), covering \ f]\ < 3.0. For 
< 1.479, the calorimeter modules are arranged in projective towers with a segmentation of 
At] X = 0.0870 X 0.0870. The ECAL and HCAL surround a tracking system with multiple 
silicon pixel and microstrip layers, covering |?/| < 2.5. Both the tracker and the calorimeters 
are immersed in a 3.8 T axial magnetic field. A detailed description of the CMS detector can be 
found inRef. [14J. 

Photons are reconstructed from clusters of energy deposited in the ECAL, using the same al- 
gorithm and granularity at the trigger level and in the offline analysis. Energy deposits within 
|A^| < 0.304 and |A?/| < 0.044 are grouped into clusters 11151 . The clustering algorithm effi- 
ciently reconstructs the energy of photons that convert in the material in front of the ECAL. 
The clustered energy is corrected taking into account interactions in the material in front of 
the ECAL and electromagnetic shower containment [16J; the correction is parametrised as a 
function of cluster size, rj, Ej, and is on average 1%. The triggers used to collect the analysed 
data sample require the presence of at least one reconstructed electromagnetic cluster with a 
minimum transverse energy of 20 or 25 GeV. The trigger is fully efficient for Ej > 21 GeV and 
\t]^\ < 1.45, defining the phase space of the measurement. Depending on the LHC instanta- 
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neous luminosity, rate-reduction factors were applied to the triggers at 20 GeV. Consequently, 
photons with Ej < 26 GeV are taken from a restricted data-set having an integrated luminosity 
of 2.1 ± 0.2 pb~^. No photon isolation criteria are applied at the trigger level. 

The event selection requires at least one reconstructed primary interaction vertex consistent 
with a pp collision [17]. The time of the ECAL signals is required to be compatible with that 
of collision products [18|. Topological selection criteria are used to suppress direct interactions 
in the ECAL APDs |19|. The residual contamination has an effect smaller than 0.2% on the 
measured cross section over the entire Ej range considered. Contamination from non-collision 
backgrounds is estimated to be negligible Iil6] . 

Photon candidates are built from ECAL energy clusters fully contained in the barrel section. 
The photon candidate pseudorapidity is corrected for the position of the primary interaction 
vertex. The absolute photon energy scale is determined with electrons from reconstructed Z- 
boson decays with an uncertainty estimated to be less than 1%. Consistent results are obtained 
with low-energy photons from decays. The linearity of the response of detector and elec- 
tronics has been measured with laser light and test beams, to a precision better than 1% in the 
energy range probed in this Letter ||T3| . Showers initiated by charged hadrons are rejected by 
requiring E^^^^/E^ < 0.05, where E^CAL sum of energy in the HCAL towers within 

R < 0.15, with = {tj — f/T)^ + {(p — (p'^Y- Electrons are rejected by requiring the absence of 
hits in the first two layers of the pixel detector that would be consistent with an electron track 
matching the observed location and energy of the photon candidate (pixel veto requirement). 

The photon candidates must satisfy three isolation requirements that reject photons produced 
in hadron decays: (1) Isotrk < 2 GeV/c, where Isojrk is the sum of the pj of tracks compatible 
with the primary event vertex in an annulus 0.04 < R < 0.40, excluding a rectangular strip 
of A// X = 0.015 X 0.400 to remove the photon's own energy if it converts into an e^e^ 
pair; (2) Isoecal < 4.2 GeV, where Isoecal is the transverse energy deposited in the ECAL 
in an annulus 0.06 < R < 0.40, excluding a rectangular strip of ^rj x = 0.04 x 0.40; and 
(3) Isohcal < 2.2 GeV, where Isohcal is the transverse energy deposited in the HCAL in an 
annulus 0.15 < R < 0.40. The requirements were designed with two other objectives in mind. 
First, the use of relatively loose photon identification and isolation selection criteria reduces the 
dependence of the results on the details of the simulation of the detector noise, the underlying 
event, and event pile-up. Second, the shape of the isolation regions is designed to allow the 
use of electrons to determine the efficiency of the isolation requirements in data. The isolation 
requirements also reduce the uncertainty on the signal due to the knowledge of the photon 
fragmentation functions. In total, 4 x 10^ photon candidates fulfill the selection criteria. 

While the isolation requirements remove the bulk of the neutral-meson background, a substan- 
tial contribution remains, mainly caused by fluctuations in the fragmentation of partons, where 
neutral mesons carry most of the energy and are isolated. A modified second moment of the 
electromagnetic energy cluster about its mean rj position, (t,^^, is used to measure the isolated 
prompt photon yield. It is calculated as 

25 25 

where Wj = max(0,4.7 + ln(£, /E)), £, is the energy of the crystal in a group of 5 x 5 centred 
on the one with the highest energy, and //; = fjj x St], where fji is the t] index of the fth crystal 
and 3ri = 0.0174; E is the total energy of the group and i] the average rj weighted by Wj in the 
same group [20J. Since cf,^^ expresses the extent in rj of the cluster, it discriminates between 
clusters belonging to isolated prompt photons, for which the cr,j,j distribution is very narrow 
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and symmetric, and photons produced in hadron decays, for which the distribution is domi- 
nated by a long tail towards higher values. Given the axial configuration of the CMS magnetic 
field, interactions with the material in front of the ECAL have a small influence on the shower 
profile along the tj direction, such that cr^rj is not affected by uncertainties on the modeling of 
such effects. The mean of the a,^,^ distributions is found to be independent of the number of 
reconstructed interaction vertices, and therefore it does not show sensitivity to pileup interac- 
tions. 

The isolated prompt photon yield is estimated with a binned extended maximum likelihood 
fit to the cr,^,^ distribution with the expected signal and background components. This is per- 
formed in each Ej bin using MiNUIT (21]. The signal component shape is obtained from pho- 
ton events generated with Pythia 6.420 p2| and the D6T parameter set [23], and simulated 
with Geant 4 [24J. The a,^,^ distribution of electrons from Z-boson decays is observed to be 
shifted when comparing data and simulated events. The shift is +(8 ± 3) x 10^^ and corre- 
sponds to 0.9% of the average of the simulated photon cr^,^ values, which are corrected for 
the observed shift. The background component shape is extracted from data by taking the 
a,j,j distribution of events in a background-enriched isolation sideband defined by requiring 
2 < Isotrk < 5 GeV/c, while keeping all other selection criteria unchanged. This choice pro- 
vides a sufficient number of events while minimising the bias to the (irjtj distribution due to 
the positive correlation between 0;^,^ and Isojrk- Both signal and background shapes are ob- 
tained separately for each Ej bin. Figure [ijillustrates the result of the two-component fit for the 
45 < £j < 50 GeV bin, which is representative of the fits in all Ej bins. The isolated prompt 
photon signal yield, N'^, is extracted with this fitting procedure. For (Tijij < 0.01, the fraction of 
isolated prompt photons in data after full selection increases from 38% at £^ = 25 GeV to 80% 
at E^ = 100 GeV. 
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Figure 1: Measured cr,jij distribution for photons with 45 < £j < 50 GeV. The fit result (solid) 
and the background component (dashed) are also shown. 



The differential cross section as a function of Ej is defined as 



d^a/dE^dr]^ = W/ {LU e AE:^ At]^), 
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where AEj is the size of the Ej bin, At]'^ = 2.9, L is the integrated luminosity, and U denotes 
bin-by-bin corrections that account for Ej reconstruction effects and finite detector resolution 
in t]^ and isolation quantities. The overall efficiency e is the product of the photon trigger, 
reconstruction, and selection efficiencies. The trigger is fully efficient for Ej > 21 GeV and 
|;/'^ I < 1.45, as previously mentioned. The efficiency of the isolation criteria is measured in data 
using an electron sample from Z-boson decays and is found to be higher than in simulation 
by pe = 1.035 ± 0.017 (stat + syst). The photon reconstruction and selection efficiencies are 
determined from Pythia events with prompt photons and are scaled by pe- The estimated 
efficiency is e = 0.916 ± 0.034 (stat + syst) and does not change appreciably with Ej or t]'^. 
Using events generated with Pythia, the values of U are calculated as a function of Ej for 
prompt photons with \f]'^\ < 1.45 and particle-level isolation less than 5 GeV. The latter is 
defined as the sum of the pj of simulated particles within R < 0.4. The resulting values of U 
decrease from 1.01 to 0.97 as Ej increases and are listed in Table [l] 

The total systematic uncertainty on the measured cross section includes contributions from the 
uncertainties in the shapes of the o;j,j distributions of signal and background, the efficiency, the 
photon energy scale, the binning of the a^,j distributions, and the U corrections. The largest 
contribution is due to the limited knowledge of the background component shape, which af- 
fects the measurement for two reasons. First, photon candidates selected from the isolation 
sideband have more associated activity in the isolation region than the true background. This 
effect is investigated by comparing the sideband and true cr,jfj distributions in simulated di-jet 
events. Events from the sideband emphasize the tail of the background cr^ri distribution, such 
that the cross section values extracted using the true background a,j,j distribution are system- 
atically lower by 15% for Ej < 85 GeV and 7% otherwise. Second, the sideband requirements 
also select some prompt photons. This effect is investigated by comparing the isolation side- 
band Crj,] distributions of simulated samples with and without prompt photons. Samples with 
prompt photons enhance the peaking part of the background distribution, such that the cross 
section values extracted using the samples without prompt photons are systematically higher 
by 12%. These two effects are checked with data by changing the isolation sideband limits so 
as to accentuate each of them. The observed variations in the extracted cross section agree with 
the estimated systematic uncertainties given above. The systematic uncertainty on the cross 
section due to the efficiencies is ±3.8%, independent of Ej and is dominated by the uncertainty 
in the efficiency of the pixel veto requirement. The full inefficiency of the pixel veto require- 
ment, estimated with simulated events, is assigned to the systematic uncertainty and is mostly 
due to the rejection of prompt photons that convert in, or before, the first layer of the pixel 
detector. The use of simulation to estimate this inefficiency is supported by the 10% accuracy 
with which the material distribution is known |25|. All the other sources of uncertainty have 
an effect on the measured cross section smaller than ±3%. 

The measured isolated prompt photon cross section as a function of Ej, including both statis- 
tical and total systematic uncertainties, is reported in Table [ij The 11% overall uncertainty on 
the integrated luminosity is considered separately. The data are shown in Fig.|2| together with 
next-to-leading order (NLO) pQCD predictions from Jetphox 1.1 [26J using the CTIO PDFs [IJ 
and the BFG set II of fragmentation functions (FF) [27J. The renormalization, factorization, and 
fragmentation scales (^r, ftp, and fif) are all set to Ej. The hadronic energy surrounding the 
photon is required to be at most 5 GeV within R < 0.4 at the parton level. To estimate the 
effect of the choice of theory scales in the predictions, the three scales are varied independently 
and simultaneously between Ej/2 and 2Ej. Retaining the largest variations the predictions 
change by (+30, —22)% to (+12, —6)% with increasing Ej. The uncertainty on the predictions 
due to the PDFs is estimated from the envelope of predictions obtained using three global-fit 
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Table 1: Isolated prompt photon cross section for \f]'^\ < 1.45 and in bins of Ej.. Uncertainties 
in the cross sections are statistical. An additional 11% luminosity uncertainty is not included 
in the systematic uncertainty (third column). The last column reports the corrections for finite 
detector resolution. A correction to account for extra activity (C = 0.97 ± 0.02) is applied to the 
theoretical predictions, as explained in the text. 
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Figure 2: Measured isolated prompt photon differential cross section and NLO pQCD predic- 
tions, as a function of Ej. The vertical error bars show the statistical uncertainties, while the 
shaded areas show the statistical and systematic uncertainties added in quadrature. A correc- 
tion to account for extra activity (C = 0.97 ± 0.02) is applied to the theoretical predictions, as 
explained in the text. The 11% luminosity uncertainty on the data is not included. 
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parametrizations, CTIO, MSTW2008 [3l, and NNPDF2.0 as recommended by the PDF4LHC 
working group [28]. This uncertainty is about ±6% over the considered Ej range. Predictions 
obtained using the CTEQ6.1M PDFs [29\, extensively used in previous comparisons with data, 
are consistent with those obtained with CTIO to within 3%. Finally using the BFG set I of FFs 
instead of the BFG set II yields negligible differences in the predictions. The theoretical predic- 
tions include an additional correction factor C{Ej) to account for the presence of contributions 
from the underlying event and parton-to-hadron fragmentation, which tend to increase the 
energy in the isolation cone. Using simulated Pythia events, C is determined as the ratio be- 
tween the isolated fraction of the total prompt photon cross section at the hadron level and the 
same fraction obtained after turning off both multiple-parton interactions and hadronization. 
Four different sets of Pythia parameters (Z2 [30], D6T, DWT, and Perugia-0 [31]) are consid- 
ered. The value C = 0.97 ± 0.02 is taken as the correction, its uncertainty covering the results 
obtained with the different Pythia parameter sets. 

As expected, the correction reduces the predicted cross section, since the presence of extra 
activity results in some photons failing the isolation requirements. 
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Figure 3: Ratio of the measured isolated prompt photon differential cross section to the NLO 
pQCD predictions. The vertical error bars show the statistical uncertainties, while the shaded 
areas show the statistical and systematic uncertainties added in quadrature. The 11% luminos- 
ity uncertainty on the data is not included. The two sets of curves show the uncertainties on the 
theoretical predictions due to their dependency on the renormalization, factorization, and frag- 
mentation scales, and on the PDFs. A correction to account for extra activity (C = 0.97 ± 0.02) 
is applied to the theoretical predictions, as explained in the text. 



Predictions from NLO pQCD are found to be in good agreement with the measured cross sec- 
tions, as shown in Figs. |2] and |3] The measured pattern is better described by the theoretical 
predictions than in previous measurements at lower ^ys and higher Xj = 2Ej/ ^Js Il8ll9l l32]]37| . 

In conclusion, a measurement of the differential cross section for the production of isolated 
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prompt photons with 21 < Ej < 300 GeV and \r]^\ < 1.45 in pp collisions at ^/s = 7 TeV has 
been presented. This measurement is performed in the kinematic regime 0.006 < xj < 0.086, 
probing a previously unexplored region at low Xt, and agrees with NLO pQCD predictions in 
the whole xj range. This measurement establishes a benchmark for photon identification and 
background estimation, and constrains the rate of one of the background processes affecting 
searches for new physics involving photons. 
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